%nonsingular part of the kernel

function [y S_m_mod] = func_Km_R_total(zeta, m, CONSTS, p_m_vec)

    K_R_nonmod        = func_Km_R_nonmod(zeta, m, CONSTS);
    
    [K_R_mod F_sum_m] = func_Km_R_mod_new(zeta, m, CONSTS, p_m_vec, true, false);

    y = K_R_nonmod + K_R_mod;
    S_m_mod = F_sum_m; % regular term in the singular part of the K_mod

end